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ABSTRACT 

The nonlinear evolution of the magnetorotational instability (MRI) in weakly 
ionized accretion disks, including the effect of the Hall term and ohmic dissipation, 
is investigated using local three-dimensional MHD simulations and various initial 
magnetic field geometries. When the magnetic Reynolds number, Reu = v a/v^ 
(where va is the Alfven speed, r] the magnetic diffusivity, and the angular frequency), 
is initially larger than a critical value R&m, crit, the MRI evolves into MHD turbulence 
in which angular momentum is transported efficiently by the Maxwell stress. If 
ReM < Re m, crit , however, ohmic dissipation suppresses the MRI, and the stress is 
reduced by several orders of magnitude. The critical value is in the range of 1 - 
30 depending on the initial field configuration. The Hall effect does not modify the 
critical magnetic Reynolds number by much, but enhances the saturation level of 
the Maxwell stress by a factor of a few. We show that the saturation level of the 
MRI is characterized by v Az /t]Q, where v\ z is the Alfven speed in the nonlinear 
regime along the vertical component of the field. The condition for turbulence and 
significant transport is given by v Az /rjQ ^ 1, and this critical value is independent of 
the strength and geometry of the magnetic field or the size of the Hall term. If the 
magnetic field strength in an accretion disk can be estimated observationally, and the 
magnetic Reynolds number v\/rjQ, is larger than about 30, this would imply the MRI 
is operating in the disk. 



Subject headings: accretion, accretion disks — diffusion — instabilities — MHD - 
turbulence 



1. INTRODUCTION 

The nonlinear regime of the magnetorotational instability (MRI) (Balbus & Hawley 1991) 
can strongly affect the structure and evolution of accretion disks. In ideal MHD, the MRI initiates 
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and sustains MHD turbulence in which angular momentum is transported outward by Maxwell 
(magnetic) stress. Thus, the MM is thought to be the most promising source of anomalous 
viscosity in disks. In weakly ionized disks, however, the coupling between the gas and magnetic 
field may be so poor that nonideal MHD effects must be considered. 

When nonthermal processes (such as irradiation by cosmic rays or high energy photons) 
dominate the ionization rate in the disk, the abundance of charged particles decreases as the 
number density of the neutral gas n n increases. At high densities (n n ^ 10 18 cm -3 ), ohmic 
dissipation dominates the evolution of the MRI (Jin 1996; Fleming, Stone, & Hawley 2000; Sano 
& Inutsuka 2001). At low densities (n n <; 10 13 cm -3 ), ambipolar diffusion dominates (Blaes & 
Balbus 1994; Hawley & Stone 1998). However, at intermediate densities, the ions are decoupled 
from the magnetic field and can drift relative to the electrons (which remain frozen- in to the field). 
Thus, in this regime Hall currents can significantly alter the MHD of the plasma, and the Hall 
term dominates the other nonideal MHD effects. Detailed calculations reveal that the Hall term 
could be important in dwarf nova disks in quiescence, and in protoplanetary disks around young 
stellar objects (Sano & Stone 2002, hereafter Paper I). 

The properties of the MRI are strongly affected by the Hall term (Wardle 1999; Balbus & 
Terquem 2001). In Hall MHD the critical wavenumber and maximum growth rate of the MRI 
both depend on the direction of the magnetic field with respect to the angular frequency vector 
fl. The Hall term can increase the maximum growth rate when the disk is threaded by a uniform 
vertical field in the same direction as i?, whereas the MRI can be completely suppressed if the 
field is oppositely directed to f2. 

In Paper I, the effect of the Hall term on the nonlinear evolution of the MRI was investigated 
using axisymmetric numerical MHD simulations. These calculations included ohmic dissipation 
as well as the Hall effect, because at some densities both processes may be important. In two 
dimensions (2D), depending on the relative amplitude of the Hall and ohmic dissipation terms 
in the induction equation, the MRI evolves into either a two-channel flow without saturation, or 
MHD turbulence that eventually dies away. 

In this paper, we continue our study of the Hall effect on the the MRI using fully three- 
dimensional (3D) numerical MHD simulations. Previous studies have shown that only in 3D 
is sustained MHD turbulence generated by the MRI (e.g., Hawley, Gammie, &; Balbus 1995). 
Moreover, the effect of nonaxisymmetric modes on the saturation amplitude and resulting stress 
can only be explored in 3D. Previous 3D simulations including only ohmic dissipation have 
shown that there exists a critical value for the magnetic Reynolds number -ReM.crit for significant 
turbulence and stress to be generated in the saturated state (Fleming et al. 2000; Sano &; Inutsuka 
2001). Moreover, this critical value depends on the field geometry in the disk (Fleming et al. 2000). 
The value of ReM,crit has important implications for the structure and evolution of accretion disks 
(Gammie 1996; Glassgold, Najita, k, Igea 1997; Sano et al. 2000). For example, the possibility 
of layered accretion and the size of the putative dead zone in protoplanetary disks depends on 
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the critical value ReM,cnt (Fleming & Stone 2002). In addition to studying the properties of the 
nonlinear regime of the MRI in 3D Hall MHD, an important goal of this paper is to investigate 
whether the inclusion of the Hall term significantly modifies the value of ReM, cr it- 

In a weakly ionized plasma, the induction equation should include terms which represent 
three non-ideal MHD effects: the Hall effect, ohmic dissipation, and ambipolar diffusion. The 
importance of these terms relative to the inductive term is determined mainly by the magnitude 
of the neutral density, the ionization fraction, and the field strength. In paper I, we calculated 
the ratios of these terms to the inductive term by solving the ionization equilibrium equations in 
both dwarf nova and protoplanetary disks. At the typical density of a dwarf nova disk (n n ~ 10 18 
cm -3 ), the ambipolar diffusion term is much smaller than the inductive term, while the Hall and 
ohmic dissipation terms are equally important and of order the inductive term. These two terms 
are also important in the inner, dense regions of protoplanetary disks. Thus, to study the MRI in 
these systems, we solve the induction equation including the Hall and ohmic dissipation terms, 
but neglecting the ambipolar diffusion term. 

Several definitions of the magnetic Reynolds number are possible; in this paper we define 
ReM = v\/i]Vi, where va is the Alfven speed, r\ is the magnetic diffusivity, and Q is the angular 
frequency. This definition uses the wavelength of the fastest growing mode of the MRI as the 
typical length scale, i.e. L = v\/Q, and therefore directly captures the effect of resistivity on the 
linear dispersion relation (Jin 1996; Sano & Miyama 1999), as well as the local properties of the 
saturated state (Paper I). Note, however, that with this definition Reu depends on the magnetic 
field strength. Moreover, it is different from the magnetic Reynolds number Re' M = Cg/rjQ, used in 
Fleming et al. (2000) by a factor v\/c 2 s , where c s is the sound speed. 

The plan of this paper is as follows. In §2, the numerical method and initial conditions used 
in the calculations are described. The results from simulations using an initially uniform vertical 
field are discussed in §3, from simulations using an initially vertical field with zero-net flux in §4, 
and from simulations using an initially uniform toroidal field in §5. The criteria for significant 
angular momentum transport, and the application of the results to protoplanetary disks, are 
discussed in §6, and our results are summarized in §7. 

2. NUMERICAL METHOD 

We adopt the local shearing box approximation (Hawley et al. 1995) for our calculations. 
The MHD equations are written in a local Cartesian frame of reference (x, y, z) corotating with 
the disk at angular frequency Q, where x is oriented in the radial direction, y is in the azimuthal 
direction, and z is in the vertical direction. Vertical gravity is ignored in this analysis, since our 
computational domain represents a region much smaller than the thickness of the disk. The gas is 
assumed to be partially ionized, and composed of ions, electrons, and neutrals. Charge neutrality 
is assumed, so that n{ = n e , where rii and n e are the number density of ions and electrons, 
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respectively. The equations we solve are then 

dp 



. )t + v-Vp = -pV v , (1) 
dv VP J x B _ - 

— + v ■ Vv = 1 2fi x v + 2qU 2 xx , (2) 

at p cp 

de „ PV • v Ann J 2 

+ v . Ve = + ; ( 3 ) 

at p &p 

dB „ / „ inn J J x B\ 

-dt =Vx ( vxB -^ ^r) = (4) 

where t» is the neutral velocity, e is the specific internal energy, 

J = ^(VxS) (5) 

is the current density, and c is the speed of light. The equation of motion (eq. 0) includes the 
Coriolis force, —217 x v, and the tidal expansion of the effective potential, 2q£l 2 x, with q = 3/2 
for a Keplerian disk. The pressure is given by P = (7 — l)pe with 7 = 5/3. The last two terms in 
the induction equation (^) represent ohmic dissipation and the Hall effect, respectively, where r\ is 
the magnetic diffusivity and e the elementary electric charge. 

These equations are solved using a finite-difference code (Sano, Inutsuka, & Miyama 1999). 
The hydrodynamics module of our scheme is based on the second order Godunov scheme (van Leer 
1979) using a nonlinear Riemann solver modified to account for the effect of tangential magnetic 
fields. The evolution of magnetic fields is calculated with the constrained transport (CT) method 
(Evans & Hawley 1988), which guarantees the divergence free condition, V ■ B = 0, is satisfied 
throughout the calculation. Each term of the electromotive force in the induction equation (||) is 
solved by an operator split solution procedure (see Paper I). 

The initial physical quantities are assumed to be spatially uniform (p = pq and P = Po) 
except for the Keplerian shear flow v y o = —qQx. The initial magnetic field is very weak for all 
models, so that radial force balance in the initial state is between Coriolis and tidal forces. Three 
initial magnetic field configurations are considered in this paper; a uniform vertical field B z = Bq, 
a zero-net flux vertical field B z (x) = BQsm(2irx/L x ), where L x is the size of the computational 
domain in the radial direction, and a uniform toroidal field B y = Bq. 

The calculations are performed in a local volume bounded by x = ±P/2, y = ±2H, and 
z = ±H/2, where H = (2/j) 1 ^ 2 c s q/Q is the scale height of the disk and c s q is the initial sound 
speed. Most of the runs use a standard grid resolution of 32 x 128 x 32 uniform zones. In the 
azimuthal and vertical directions, periodic boundary conditions are used. For the radial boundary, 
a sheared periodic boundary condition (Hawley & Balbus 1992) is adopted. The magnetic flux 
within the shearing box must be conserved unless a net flux of radial field exists (Hawley et al. 
1995). However, we have found that numerical errors in the net flux of the vertical and azimuthal 
fields are non-negligible (more than 10 % in 10 orbits) when the EMF in the ghost zones at the 
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radial boundary are constructed from applying the sheared periodic boundary condition to the 
magnetic field and velocity. Instead, if we apply the the sheared periodic boundary condition to 
the azimuthal component of the EMF, the error in the vertical flux is reduced to ~ 0.1 %. 

We choose the normalizations po = 1, H = 1, and Q = 10~ 3 . Then the sound speed and gas 
pressure are initially 2c 2 /7 = 10~ 6 and Pq = 5 x 10~ 7 , respectively. Initial perturbations are 
introduced as spatially uncorrelated pressure and velocity fluctuations. These fluctuations have a 
zero mean value with a maximum amplitude of \6P\/P = W~ 2 and \Sv\/c s0 = 10" 2 . 

Our calculations are characterized by three dimensionless parameters. The first is the ratio of 
the gas to initial magnetic pressure, 



8vrP _ 2c 2 
B o ~ l v \o 



Po = ^ = , (6) 



where vao = Bq/(Atipq) 1 / 2 is the initial Alfven speed. This parameter measures the initial field 
strength. The second is the magnetic Reynolds number, which is defined as 

fle„o = f . (7) 

We assume the magnetic diffusivity r\ is constant in our calculations. This parameter measures 
the importance of ohmic dissipation. Finally, the third parameter is 

_ C B n 

2iren e0 v A0 

where n e Q is the initial number density of electrons. We assume the electron abundance is constant 
throughout our calculations, thus n e = n e Qp/pQ. The value of Xq measures the importance of the 
Hall term; note that the sign of Xq depends on the direction of the magnetic field. 



3. SIMULATIONS WITH A UNIFORM VERTICAL FIELD 

Simulations that begin with a uniform vertical field are listed in Tables [l] - 13[ The critical 
wavelength A cr it and the maximum growth rate o" max are obtained from the dispersion relation 
derived by Balbus & Terquem (2001) for axisymmetric perturbations. Disturbances with a 
wavelength longer than A cr it are unstable to the MRI. The critical wavelength is of the order 
of 2ttva/Q in most cases when ohmic dissipation is inefficient (Re mo > !)• However, when 
the Hall parameter is just below zero (—4 < Xq < 0), A cr i t — > and small-scale perturbations 
become unstable to the MRI. Ohmic dissipation, on the other hand, can stabilize small-scale 
perturbations, and thus make A cr it longer and cr max smaller when ReMO < 1. In 2D, when Xq < 
or ReMO < 1 saturation of the instability occurs, whereas exponential growth of a two-channel flow 
without saturation persists for other models (Paper I). In 3D without the Hall term, the nonlinear 
evolution of the MRI is characterized by recurrent appearance of two-channel flow (Fleming et al. 
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2000; Sano & Inutsuka 2001). We now examine how this evolution is changed by the inclusion of 
the Hall effect. 

In addition to model parameters, Tables |l| — 13| contain several fundamental quantities that are 
referred to in the following discussion. Hereafter, the single brackets (/) imply a volume average 
of quantity /, whereas double brackets ((/)) denotes a time and volume average. If not otherwise 
stated, the time average is taken over the last 10 orbits of the calculation for the models in this 
section. 



3.1. The Fiducial Models 

We describe three models in detail: models Z2 (Xq = 2), Z3 (Xq = 0), and Z4 (Xq = —2). 
These fiducial models have the same initial field strength and magnetic Reynolds number, 
Po = 3200 and Rcmo = 100. The only difference between the non-zero Xq models (Z2 and Z4), is 
the direction of the vertical field. Ohmic dissipation is so weak in these models that the evolution 
is dominated by the Hall effect. 

Figure [j] shows the time evolution of the volume-averaged magnetic energy (B 2 /8tt)/Pq for 
these models, where time is measured in orbits i ro t = 2ir/Vt. During the linear phase of the MRI, 
magnetic energy is amplified exponentially. As expected from the linear analysis, the Xq = 2 
model has a larger growth rate than the Xq = —2 run initially. In the nonlinear regime, MHD 
turbulence is sustained in all three models. The saturated magnetic energy shows large fluctuations 
whose amplitude depends on the size of the Hall parameter: as Xq decreases, the amplitude of the 
variability becomes smaller. The saturation level of the magnetic energy is higher in the positive 
Xq model than in the negative Xq model. 

The efficiency of angular momentum transport is given by the turbulent stress, 

w xy = + P°x° v v > ( 9 ) 

where the first and second terms are the Maxwell stress (wm = —B x B y /4nr) and Reynolds stress 
(wr = pv x 5v y ), respectively. The total stress w xy is related to the a parameter of Shakura & 
Sunyaev (1973) by a = w xv /Pq. Figure |2| shows the time evolution of the Maxwell stress normalized 
by the initial pressure Pq for the fiducial models. The time-averaged Maxwell and Reynolds 
stresses are listed in Table |l]. In MHD turbulence generated by the MRI, the Maxwell stress is a 
few times larger than the Reynolds stress. Large time variability can be seen in the evolution of 
the stress, with the amplitude of fluctuations in the Xq = 2 run about A(wm) /Pq ~ 0.4, which 
is much larger than the time-averaged value of the stress ({wm))/Po = 0.16. In the model with 
negative Xq, on the other hand, the variation is much smaller throughout the evolution. 

The spike-shaped variations in the magnetic energy and stress correspond to the recurrent 
appearance and breakup of the two-channel flow. These same variations could be seen in the 
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3D simulations of the MRI without the Hall effect studied by Sano & Inutsuka (2001). Figure || 
shows images of the magnetic energy distribution in the Xq = 2 run (Z2) at 18, 21, and 25 orbits. 
The velocity field is also shown by arrows. The top and bottom panel are chosen at times near 
the peak of spikes in the magnetic energy, whereas the middle panel is near a minimum. The 
large-scale structure of the flow in the top and bottom panels clearly shows the axisymmetric 
two-channel flow. The spatial dispersion in the magnetic pressure is very high during these 
phases; at 18 orbits (£P mag } 1/2 / '(-Pmag) = 1.5. The magnetic pressure is comparable to the gas 
pressure ((-P m ag)/(-P) = 0.65), and the gas pressure and density have large spatial dispersions 
{{8P 2 ) X ' 2 /{P) = 0.52 and (fy 2 ) 1 / 2 /^) = 0.38). 

The axisymmetric channel flow is unstable to nonaxisymmetric modes of the parasitic 
instability (Goodman & Xu 1994). After the breakup of the channel flow, organized large-scale 
structure disappears and the disk is occupied by disorganized MHD turbulence (Fig. [3], middle 
panel). The magnetic energy at this phase is an order of magnitude smaller than at the peak of 
the spikes, so that (-P)/(-P m ag) = 15 at 21 orbits. The magnetic pressure is still spatially highly 
fluctuating, {SP^g) 1 ^ 2 / (P ma ,g) = 1.1, but this fluctuation is small relative to the gas pressure, 
(5P mag ) 1 / 2 /(P) = 0.075. Since the density fluctuation decreases as the fluctuation in the magnetic 
pressure relative to the gas pressure (8P^ ag ) 1/2 /(P) decreases (Turner, Stone, & Sano 2002), the 
spatial dispersion in the gas pressure and density is small at this phase; (5P 2 ) 1 / 2 / (P) = 0.14 and 
(8p 2 )V 2 /{p) = 0.085. 

The snapshots of the magnetic energy in the negative Xq model (Z4) is shown in Figure ||. 
In contrast to the Xq = 2 run (Z2) shown in Figure ||, this model has little time variability. To 
allow direct comparison, the contour levels are the same in both figures. The magnetic energy is 
typically comparable to or smaller than that during low energy phase of the Xq = 2 run. The 
emergence of the channel flow cannot be seen, and disorganized MHD turbulence is sustained 
throughout the evolution. Because the critical wavenumber fc cr it — > oo in this case, small-scale 
disturbances are noticeable in all directions. The gas pressure is always much larger than the 
magnetic pressure and the ratio of their time average is ((-P)}/({P m ag)) = 32. The spatial dispersion 
in the magnetic energy is very large, {SP^g) 1 ^ 2 / (-Pmag) = 1-5, but the gas pressure and density are 
almost uniformly distributed, {5P 2 ) 1/2 /{P} = 0.093 and {5p 2 ) 1/2 /{p} = 0.054. These quantities 
are quite similar to those during the low magnetic energy phases of the Xq = 2 run. 

To compare the characteristics of MHD turbulence in the fiducial models, Table § lists a 
number of time- and volume-averaged quantities. Although the magnetic and kinetic energy in 
the positive Xq run are about 4 times larger than those in the negative Xq run, interestingly the 
turbulence of three fiducial runs has many similarities. For example, the magnetic and perturbed 
kinetic energy are almost equal for all the models, with the ratio {{p5v 2 / '2)) / {(B 2 /8ir)) ~ 0.4. 
Because of the shear motion, the azimuthal component of the magnetic field is amplified most 
efficiently. Each component of the magnetic field energy has a similar ratio for all cases, 
((B 2 )) : {{B 2 )) : {{B 2 )) f» 3 : 20 : 1. The perturbed kinetic energy is slightly anisotropic, 
((pv 2 J : ((p5v 2 )) : ((pv 2 )} » 3 : 2 : 1. 
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The model parameters, /3q, ReMo, and Xq, are denned by the initial state, but they are 
variable in time. The time-averaged values of these parameters in the fiducial models are listed 
in Table |4j. The magnetic pressure is smaller than the gas pressure even in the nonlinear stage, 
so that the time- and volume-averag ed plasma beta ((/?)) = {{P/(B 2 /8tt))) is of the order of 100. 
However, the gas pressure is increasing linearly throughout the evolution due to the dissipation 
of magnetic field. Then, the plasma beta {{(3)) must be increasing in time unless some cooling 
processes are included. When the magnetic field is amplified, the magnetic Reynolds number 
becomes larger than the initial value. The time-averaged value is ((Rcm)) = ((^a/ 7 ?^)) > 10 4 f° r 
the fiducial models, so that the ohmic dissipation at the nonlinear stage is less efficient. We define 
the effective Hall parameter in the nonlinear regime as 

x ca = cBzQ = 2cpBzQ (10) 

c 2iren e v\ en e B 2 

Because the Hall parameter is inversely proportional to the field strength, the volume-averaged 
value K-Xeff)! decreases as the field is amplified. The time-averaged value is {{X e g)) = 0.003 and 
—0.03 for the Xq = 2 and —2 run, respectively. Therefore, the effect of the Hall term may be 
reduced in the nonlinear regime for these cases. 



3.2. Saturation Level 

In this subsection, we explore in more detail the influence of the model parameters on the 
time- and volume-averaged stress in the saturated state. Time-averaged Maxwell and Reynolds 
stress are listed in Table [l] - 13| for all models. As shown by the fiducial models, the stress exhibits 
large time variability in the turbulent state. For all the models, the time averaging is taken over 
the last 10 orbits of the calculation, which is longer than the typical timescale of variations in the 
stress. The time-averaged quantities taken over the last 40 orbits of the fiducial models (Z2 and 
Z4) are similar to the 10 orbit averages, and the difference in the Maxwell and Reynolds stress is 
less than 10 %. The time-averaged magnetic Reynolds number defined by the vertical magnetic 
field ({v\ z /r]Q,}} is also listed in the tables. The meaning of ((v\ z /r]Q)) is discussed later in §6.1. 



3.2.1. Effect of Initial Field Strength 

First we study the effect of the initial field strength. Figure [| depicts the saturated Maxwell 
stress as a function of the Hall parameter Xq for various models with different [3q. All the 
models shown in this figure has the same magnetic Reynolds number Re mo = 1- For any (3q, 
the saturation levels in the positive Xq runs are higher than those in the negative Xq runs. The 
ratio of the stress between the Xq = ±2 runs is 28, 4.1, and 2.6 for (3 = 800, 3200, and 12800, 
respectively. All models with Xq > show large time variability due to the nonlinear growth of 
the channel flow, while no growth of the two-channel flow can be seen in all models with Xq < 0. 
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If we compare models with the same Xq, larger (3q models have a lower saturation level, which 
means that the magnetic energy and stress increase as the initial field strength increases, as in the 
ideal MHD cases (Hawley et al. 1995). 

Since the linear growth rate of the MRI is higher for Xq > (Balbus & Terquem 2000; Paper 
I), this may account for the higher saturation level in this case. In addition, the evolution of the 
MRI shows the recurrent growth of the channel flow when the Hall parameter is Xq > 1. Since the 
channel flow can amplify the magnetic field more efficiently than disorganized MHD turbulence, 
this could also be a reason for the larger saturation level in the positive Xq runs. This result may 
also be understood in terms of the linear properties of the MRI: If Xq > 0, the critical wavelength 
is proportional to the field strength A cr it ~ va/&, so that A cr it increases as the magnetic energy 
is amplified, leading to the emergence of large-scale channel flows. When the Hall parameter is 
negative, on the other hand, the critical wavenumber for the MRI is infinity, so that small-scale 
perturbations are unstable. The MRI continuously excites small-scale disturbances in this case, 
and these small fluctuations impede the nonlinear growth of the two-channel flow. 

3.2.2. Effect of Magnetic Reynolds Number 

When the magnetic Reynolds number is very small, ohmic dissipation can dramatically reduce 
the linear growth rate (Jin 1996) and the nonlinear saturation level of the MRI (Sano & Inutsuka 
2001). The dependence of the saturated stress on the magnetic Reynolds number is illustrated in 
Figure ||, which shows the time-averaged stress for the models with ReMo = 100, 1, and 0.1. The 
same field strength is used for all the models in this figure ((3q = 3200). We find that the positive 
Xq runs always have a larger stress than the negative Xq runs. The differences are by a factor 
of 4 - 6 in the ReMo = 100 and 1 runs. For very resistive models with Remo = 0.1, the ratio of 
the stress between the Xq = 4 and —2 runs is about 100. However the stress is of the order of 
10~ 5 - 10 -4 , and this is more than 2 orders of magnitude smaller than the less resistive cases 
(ReMQ > 1). This suggests that strong dissipation can weaken the turbulence even when the Hall 
term is included. 

3.2.3. Effect of Resolution and Box Size 

Table |2| lists the model parameters and saturation levels for high resolution runs designed 
to study the effect of numerical resolution. Except for the grid resolution, the model parameters 
of Z2H and Z4H are the same as those of the non-zero Xq fiducial models Z2 and Z4. The high 
resolution models are followed only to 10 orbits, so that the time average is taken over the last 5 
orbits. In the standard resolution models (Z2 and Z4), the time-averaged stress over every 5 orbits 
after 10 orbits has up to 70 % difference compared to the average through the last 40 orbits. It 
may be that the time-averaged stress in the high resolution models has uncertainty of the similar 
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size. We find that most of the nonlinear features shown by the fiducial models are independent 
of the grid resolution. The magnetic energy and stress in the Xq = 2 run (Z2H) is larger than 
those in the Xq = —2 run (Z4H) by a factor of about 4. The saturation levels in models Z2H and 
Z4H are close to those in the standard resolution cases. The positive Xo run shows the nonlinear 
growth of the channel flow, and the negative Xo run evolves into disorganized MHD turbulence 
with many small-scale fluctuations. 

Table ^ lists the models computed with a larger computational box than the fiducial models, 
namely 2H x 8H x 2H. The grid resolution is 64 x 256 x 64, so that the grid spacing is the same 
as the standard models listed in Table |l|. The other parameters of models Z2L and Z4L are the 
same as the fiducial models Z2 and Z4. As in the standard box size cases, the magnetic energy and 
stress in the positive Xq run (Z2L) is larger than those in the negative Xq run (Z4L) throughout 
the evolution. The frequency of spike-shaped variations in the stress is reduced with a larger box 
size, and thus disorganized MHD turbulence lasts most of the time. The stress in the large box 
model with Xo = 2 (Z2L) is comparable to that in the standard model Z2, but the Xo = —2 model 
(Z4L) is twice as large as the standard model Z4. Although the saturation level for models with a 
uniform vertical field may be proportional to the box size for ideal MHD simulations (Hawley et 
al. 1995), our sample of simulations is too small to confirm this dependence in this study. 



3.3. Characteristics of Saturated MHD Turbulence 

It is of interest to study the properties of the MHD turbulence driven by the MRI in Hall 
MHD. Figure 0a and 0b show the Fourier power spectra of the magnetic energy along the k x , k y , 
and k z axes at 25 orbits for the fiducial models Z2 (Xo = 2) and Z4 (Xo = —2), respectively. 
The spectra are averaged over 10 snapshots within 0.1 orbits. In the Xo = 2 run, the large-scale 
channel flow is growing at that time. Thus, smaller wavenumbers have larger power, especially in 
k z , and the power at the inertial range declines as k~ 4 , similar to the ideal MHD case (Hawley et 
al. 1995). For the Xo = —2 run (Fig. |7|b), the amplitude of the power is smaller than the Xo = 2 
run, although the shape of the spectra is quite similar to the positive Xo run. The slope of the 
power is slightly gentler in the Xo = —2 run, probably because modes with larger k are unstable 
to the MRI in this case. The power spectra at low magnetic energy phase of the X = 2 run (which 
is dominated by disorganized MHD turbulence) are similar to those of the negative Xo model. 

When ohmic dissipation is inefficient Re mo > 1, the MHD turbulence generated by the 
MRI appears to have the characteristic properties. For example, the Maxwell stress is always 
proportional to the magnetic pressure, 

(W = (0.455 ± 0.023) ((P mag )) , (11) 

where the average and the standard deviation are taken from the 17 models with Re mo > 1 listed 
in Table 1. Note that this average includes models both with and without the Hall effect. The 
deviation of this ratio is extremely small, and this is independent on Xq. The relation given by 
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equation (11) is the same as in the ideal MHD cases (Hawley et al. 1995). For the models with 
Rcmo > 1, the ratios between each component of the magnetic and perturbed kinetic energy 
({(B 2 )) : ({By)) : ((B%)) and {(pv 2 )) : ((pSVy)) : {(pv 2 ))) take quite similar values to the fiducial 
models. These ratios have also no correlation with Xq. The average ratio of the Maxwell and 
Reynolds stress is 

{(w M )) = (AA9±im)((w R )) . (12) 

The ratio ((wm)) / {(wr)) is slightly dependent of Xq, and the average is 5.77 and 3.03 for the 
Xq = 2 and —2 runs, respectively. 

When ReMO < 1, the turbulence in the nonlinear regime is weakened by ohmic dissipation. 
In this case, the power spectra of the turbulence is quite different compared to the ReMO > 1 
runs. Figure || shows the power spectrum for model Z13 with (3q = 3200, ReMO = 0.1, and Xq = 0. 
The saturation level of the Maxwell stress is very low for this model ({(wm)) / Po = 3.1 x 10 -4 ). 
Other models with ReMO = 0.1 including the Hall term (models Zll, Z12, and Z14) show little 
difference in the energy spectra compared to the Xq = run. The typical diffusion scale defined 
as fcdiff = v/ v A0i as well as the critical wavenumber for the MRI, A; cr it, is shown on the plot of 
the power spectra of the magnetic energy (Fig. [8|, left panel). The diffusion length is close to the 
critical wavelength in this case. For scales smaller than the diffusion length, fluctuations in the 
magnetic field dissipate faster than the Alfven timescale. Therefore, the spectrum shows a very 
steep decline and the slope is proportional to k~ s in the dissipation regime k ^ fcdiff- The right 
panel of Figure || shows the power spectra of the perturbed kinetic energy, which also is a steeply 
decreasing function of k. 

We find that the energy distribution of the turbulence in the Re mo = 0.1 runs is also quite 
different from that in the less resistive models. In Table |], time- and volume-averaged quantities 
for the ReMO = 0.1 runs (Z12, Z13, and Z14) are listed. The saturated level of the magnetic energy 
is much lower than the ReMO = 100 runs. The magnetic field is amplified by the MRI during 
the linear phase, but dies away due to the ohmic dissipation, so that in the nonlinear regime 
the magnetic energy returns to its initial value. Because the net flux through the computational 
volume is conserved, the magnetic energy can never completely die away. The energy in the 
perturbed velocity is larger than in the magnetic field, ((p5v 2 /2)) 3> ((B 2 /8ir)), and the Reynolds 
stress is a few times larger than the Maxwell stress for the Re mo = 0.1 models. The large kinetic 
energy is a remnant of the linear growth of the MRI. Although the perturbed magnetic field 
generated by the growth of the MRI is efficiently dissipated in these models by the large resistivity, 
the perturbed kinetic energy remains large due to the small viscosity. Therefore, if the disk is 
linearly unstable and the magnetic Reynolds number is small (ReMO < 1) 5 the Reynolds stress 
could dominate the Maxwell stress but the efficiency of angular momentum transport is very small 
(a ~ 10" 5 - 10" 3 ). 
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3.4. Critical Magnetic Reynolds Number 

In this subsection, we consider the critical value of the magnetic Reynolds number for 
significant angular momentum transport in accretion disks with uniform vertical fields (as shown 
by previous studies, this critical value depends on the initial field geometry, Fleming et al. 2000). 
The dependence of the saturation level on the magnetic Reynolds number ReMo is shown by 
Figure ||. Open circles and triangles denotes the results without the Hall term (Xq = 0) for the 
model with (3q = 3200 and 12800, respectively. When ReMo > 1, the saturation level is of the 
order of 0.01 and almost independent of the magnetic Reynolds number. As shown by Figure ||, 
the stress is larger when the initial field is stronger. If the magnetic Reynolds number is less than 
unity, on the other hand, the saturation level is reduced by a large factor. For the ReMo = 0.1 
runs, the stress is about 10 for both cases. Thus, the criterion for significant turbulence is 
Rcmo 1 independent of the field strength. Note that the critical value of the magnetic diffusivity 
77 depends on the field strength, because the magnetic Reynolds number is a function of the Alfven 
speed (r/ cri t = v\ /QRe M o,crit ~ v\ Q /Vt). 

Next the effect of the Hall term on the critical magnetic Reynolds number is considered. 
Filled circles and triangles show the saturated stress in the models with non-zero Hall parameter 
Xq = 4 and —2, respectively. Although the stress in the Xq = 4 runs is always several times larger 
than that in the Xq = —2 runs, the dependence on the magnetic Reynolds number is the same as 
in the models without the Hall term. That is, the saturation level is almost independent of Re\jo 
when ReMo > 1, but drops significantly as the magnetic Reynolds number decreases if ReMo < 1- 

In 2D simulations, all the models with Re mo < 1 undergo rapid decay of the magnetic energy 
both with and without the Hall term (Paper I). Thus, the dependence of the saturation level on 
the ohmic dissipation and the Hall effect is consistent with the 2D results. When ReMo > 1> the 
magnetic Reynolds number ((vj^/rjQ,)) at the nonlinear stage is larger than the initial value, because 
the magnetic field is amplified by the MRI. If ReMo < 1, the magnetic energy is unchanged by the 
instability, and thus ((v\/t]Q,)} remains less than unity. Therefore, the critical value ((v\/t]Q,)} ~ 1 
is valid even in the nonlinear regime for uniform B z models. 

4. SIMULATIONS WITH A ZERO-NET FLUX VERTICAL FIELD 

We next consider simulations that begin with a zero-net flux vertical field, that is 
B z {x) = BQs'm(2irx/L x ) where Bq is a positive constant. For this case, the Hall parameter is given 
by X(x) = Xq/ sm.(2-Kx/L x ), where Xq = 2cpo$l/ 'ere e o-E>o- Thus, the region x < has positive X 
while the region x > has negative X, and the minimum of the absolute value |A(x)| is Xq. 
Table || lists the models computed with a zero- net flux B z . The initial plasma beta (3q, the critical 
wavelength A cr it, and the maximum growth rate <r max in this table are given for B z = Bq and 
X = Xq. The saturation level of the Maxwell and Reynolds stresses and the magnetic Reynolds 
number are also listed in the table. No data in the columns for the saturation level means that it 
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is less than 10 8 . The time average is taken over the last 20 orbits for all models in this section. 



Figure 1C shows the time evolution of the magnetic stress for models SI (Xq = 0), S2 
(Xq = 2), and S3 (Xq = 4). All the other parameters for these models are identical ((3q = 3200 
and ReMo = 100). Because the magnetic Reynolds number for these models is very large, the 
Hall effect dominates the evolution of the MRI. The stress in the Xq = 4 run has large amplitude 
time variations and a higher saturation level than that in the Xq = run. Thus, the Hall effect 
enhances the saturation level of the stress, as in the uniform B z cases. The difference in the time 
average ((wm)) is a factor of 3 between the Xq = and 4 runs. This saturation level is an order of 
magnitude lower than the uniform B z models with the same (5q and ReMo- 

During the linear phase, the MRI grows most rapidly in the half region with positive X, 
because most of the region with X < is linearly stable for the MRI or has a smaller growth rate 
(see Paper I). The amplified magnetic field gradually affects the structure of the other (X < 0) 



half region, and after several orbits the entire region becomes turbulent. Figure 11 shows images 
of the magnetic energy during the turbulent phase for model S3 (Xq = 4) at 35, 40, 45, and 50 
orbits. The magnetic field vectors in x-z plane are also shown by arrows. A large time variation 



in the stress occurs near 35 orbits (see Fig. |T0| ), however from figure 11 no channel flow is evident 
at this time. Instead disorganized MHD turbulence is sustained throughout the evolution. The 
amplitude of time variation in the stress is slightly smaller compared with uniform B z cases. At 
the turbulent phase, the initial distribution of the vertical field disappears completely, and the 
positive and negative X regions are randomly distributed. The effective Hall parameter at the 
nonlinear stage is very small and negative for this case, ((X e fj)) = —0.018. At 50 orbits, the spatial 
dispersion is large in the magnetic pressure 

(( 5 - P mag) 1/2 /(- p ma g ) = 1-0) but small in the density 
and gas pressure ((6P 2 ) 1 / 2 1 '(P) = 0.12 and (<5 / o 2 ) 1/ V(P> = 0.069). These numbers are similar to 
the uniform B z cases during phases of low magnetic energy. 

Since there is no net-flux for this field configuration, the magnetic field within the shearing box 
can completely die away. In fact, the magnetic energy during the nonlinear regime is decreasing 



in time for models with small Re mo- Figure 12 shows the time evolution of the magnetic energy 
for models with Re mo = 1- Ohmic dissipation affects the linear characters of the MRI in these 
cases. When the Hall parameter is small, models S8 (Xq = 0) and S9 (Xq = 2), the magnetic 
energy is amplified by the MRI during the linear phase, but this amplified field is not sustained. 
After a few tens of orbits, the magnetic energy and stress decrease due to ohmic dissipation until 
the end of the calculation. In the Xq = 4 run (S10), however, the amplified magnetic energy is 
sustained for at least 50 orbits even with ReMo = 1- The time evolution of run S10 is quite similar 



to the less resistive models shown in Figure 10. The saturation level of the Maxwell stress in 
model S10 is {{wm))/Po = 0.033 and this is comparable to those in models S3 (ReMo = 100) and 
S6 (Re mo = 10). This enhancement is caused by a faster linear growth rate for the MRI at larger 
Xq (Paper I). Once the field is amplified, the efficiency of the Hall term and ohmic dissipation is 
reduced, and thus turbulence can be sustained. 
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Table || lists the time- and volume-averaged quantities for the models with ReMo = 100 (SI 
and S3) and R&mq = 1 (S8 and S10). Models SI and S8 do not include the Hall effect (Xo = 0), 
whereas models S3 and S10 have X = 4. The properties of the turbulence in the X = 4 runs (S3 
and S10) are nearly identical. The saturation level of the magnetic energy is ((B 2 / 8n)) / Po ~ 0.1, 
which is 3 times larger than the perturbed kinetic energy. The ratio of the Maxwell and Reynolds 
stress is ((wm)} / {(wr)} ~ 4. The ratio of each component of the magnetic energy in the turbulence 
is ((B 2 )) : {(By)) : ((B 2 )) ~ 2 : 18 : 1, and this is quite similar to that in the uniform B z models. 
Since the magnetic Reynolds number in the nonlinear regime ((ReM)) is large, ohmic dissipation is 
ineffective. Moreover, the Hall effect may also be unimportant in the nonlinear regime, because 
the effective Hall parameter is small ((X e g)) ~ —0.02. The saturated quantities in model SI are 
also close to the Xq = 4 runs. For these three models, the saturated level of a is of the order of 
0.01. 

The turbulence in model S8, on the other hand, is strongly affected by ohmic dissipation. The 
stress a = 7.8 x 10 -6 at 50 orbits is much smaller than the other runs, and is decreasing in time. 
The properties of the turbulence, such as {(B 2 )) / ((B 2 )) , are also very different from the other runs. 
The magnetic Reynolds number in the nonlinear regime is still small {(Rcm)) = 1-2, which means 
the ohmic dissipation is important throughout the evolution of this model. 



4.1. Critical Magnetic Reynolds Number 

Here we estimate the critical value of the magnetic Reynolds number ReM0,crit for zero-net 
flux B z models. Figure [l3| shows the saturation level of the Maxwell stress for the models without 
the Hall term (Xq = 0). Three different cases with (3q = 800, 3200, and 12800 are shown in 
this figure. When ReMo = 100, the stress for all the models is more than 0.01, so that angular 
momentum transport is efficient. If the magnetic Reynolds number is below a critical value, the 
stress drops dramatically. For the /?o = 3200 runs, the stress is {{wm))/Po = 0.018 and 2.0 x 10 -5 
at ReMo = 10 and 3, respectively. Thus the stress decreases about 3 orders of magnitude with a 
small difference in Re mo- 

The critical value is ReMo,crit ~ 10 for (3q = 3200, but this value is found to depend on the 
initial field strength. The condition for ((wm))/Pq > 0.01, for example, is ReMo > 30, 10, and 3 
for (3q = 800, 3200, and 12800, respectively. When ReMo > -ReMo.crit, the saturation level of the 
Maxwell stress is nearly independent of (3$ (Hawley, Gammie, &: Balbus 1996) and this is different 
from the uniform B z case. In order to initiate the turbulence, the MRI must grow faster than the 
dissipation rate of the initial field. Thus, the critical value for the instability should be given by 

A "" ReMO\/Po (t~\ ~ const. , (13) 



(where L is the length over which the initial field varies) so that the critical value of Rcmo 
is proportional to f3 1 ^ 2 L _1 . This idea is roughly consistent with the results shown in Figure 
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p~3| . If L decreases, the critical value should increase. If fact, when we use a field distribution 
B z {x) = Bq sm(6irx/ L x ) (i.e., L = L x /3), the criterion becomes Reuo > 30 for (3q = 3200. This 
suggests that if the zero-net flux magnetic field has small-scale structure initially, then the ohmic 
dissipation can suppress the initiation of the MRI with a small magnetic diffusivity r\. 

Next we examine the effect of the Hall term on the critical magnetic Reynolds number. Figure 



14 shows the saturated magnetic energy as a function of the initial magnetic Reynolds number 
ReMO- F° r comparison, the magnetic Reynolds number Re' M0 used in Fleming et al. (2000) is also 
shown in the figure. The initial field strength for all the models in this figure is [3q = 3200. Open 
circles denotes the models without the Hall term. For Xq = 0, as discussed above, the saturation 
level is almost constant for ReMO > 10, but decreases as the magnetic Reynolds number decreases 
if Re^io < 10- When Re^jo = 0.3, the magnetic energy does not show any increase during the 
evolution, and almost dies out: ((B 2 / 8tt}) / Pq ~ 10~ 12 at the end of the calculation. The critical 
value for significant turbulence, {(B 2 / 8tt)) / Pq > 0.01, is ReM,ait ~ 10. This value corresponds to 
i?e' MOci . it ~3 x 10 4 , which is consistent with the results of Fleming et al. (2000). 

For the Xq = 2 runs depicted by filled triangles, the effect of the Hall term is not large; the 
critical magnetic Reynolds number is still ReM0,crit ~ 10. However, if the size of the Hall term 
increases, we find that ReM,crit shifts to smaller values. The Xq = 4 runs are shown by filled circles 
in this figure. The saturation level in the Rej^o = 1 run is increased to ((B 2 / 8tt)} / Pq = 0.088. But 
the ReMO = 0.3 run shows no growth of the MRI. Thus, the critical value for {(B 2 /8ir}) / Pq > 0.01 
becomes -ReA/o,crit ~ 1, which is an order of magnitude smaller than the case without the Hall 
term. 

The Hall parameter Xq can take a large range of values when ReMO is small, because the 
maximum value is Xq ~ 100/i?eA/o m the Hall regime (Paper I). The negative X region (x > 0) is 
linearly stable when Xq > 4, however this is potentially offset by the larger linear growth rate in 
the regions with a large positive Hall parameter. We have calculated models with Xq = 100 (filled 
squares in Fig. 14) and 1000 (cross). For ReMO = 0.3, the saturation level is enhanced by many 



orders of magnitude when Xq = 100, but the stress is still very small ((wm))/Po = 7.0 x 10~ 5 
compared with the less resistive models (ReMO > 10). For more resistive models (ReMO = 0.1), 
the MRI does not operate even with Hall parameters as large as Xq = 100 and 1000. Note that in 
the model with ReMO = 0.1 and Xq = 1000 (S16), the critical wavelength is longer than the scale 
height of the disk, meaning such large values of Xq cannot result in enhanced transport in real 
disks. Therefore, the change in the critical magnetic Reynolds number due to the Hall term is at 
most an order of magnitude for zero-net flux vertical fields. 



5. SIMULATIONS WITH A UNIFORM TOROIDAL FIELD 

Lastly we investigate the evolution of the MRI starting with a uniform toroidal field, B y = Bq. 
Table ^ lists the models calculated with this field geometry. Axisymmetric perturbations are stable 
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if the field is purely toroidal. Since no linear dispersion relation has been obtained for the MRI for 
nonaxisymmetric perturbations including nonideal MHD effects, the characteristic length of the 
MRI defined as Amri = 2t:vaq/VL is listed in the table instead of A cr i t and cr max . The length Amri 
corresponds to the wavelength in the azimuthal direction of the most unstable mode. The table 
gives the saturation level of the Maxwell stress {(wm))/Po, the Reynolds stress ((wr))/Po, and the 
magnetic Reynolds number ((v Az /rjQ)}, if they are more than 10~ 8 . In some cases the uniform B y 
models must be evolved for very long times to reach a saturated nonlinear regime. Time averages 
are taken over the last 40 orbits for all the models in this section. 

We still use the Hall parameter Xq defined by equation (||) as a model parameter. Balbus & 
Terquem (2001) examined the Hall effect on the nonaxisymmetric behavior of linear perturbations. 
The critical wavenumber for the MRI becomes longer or shorter depending on the sign and size of 
the Hall parameter. The effective Hall parameter, Ha, for nonaxisymmetric disturbances with a 
wavenumber k is given by 

H _ c(fc • B)(k • O) _ ^ / k y v Ao \ /&zUao\ 

Thus, the parameter Ha depends on the wavenumber of each mode. We begin the simulations 
with small random perturbations, so that modes with every allowed k are included initially. We 
expect the most unstable mode to dominate the linear phase of the MRI. In the ideal MHD 
limit for a uniform B y , the most unstable mode has a characteristic wavenumber k y ~ vao/Q in 
the azimuthal direction, but no scale in the vertical direction k z — > oo (Balbus & Hawley 1991). 
In the numerical calculations, the grid resolution constrains the minimum length in the vertical 
direction, which is typically k z VAo/Q ~ 10. Because the parameter Ha is larger than Xq by the 
factor k z VAo/Q, we expect the Hall term will have an effect even with small Xq. 

Figure [IB] shows the time evolution of the Maxwell stress for the models with Xq = 0, 0.2, and 
0.4, which are models Y2, Y3, and Y4, respectively. All the models have the same field strength 
and the same magnetic Reynolds number, [5q = 100 and ReMO = 100. The evolution is dominated 
by the Hall effect because ohmic dissipation is inefficient (ReMO = 100 3> 1). We find the timescale 
to reach saturation is very sensitive to the Hall parameter. Without the Hall term (Xq = 0) it 
requires 15 orbits, but in the Xq = 0.4 run the turbulent state begins at about 8 orbits. Thus, as 
discussed above, the linear phase of the nonaxisymmetric MRI is affected by even small Xq. 

Despite the sensitivity of the linear growth rates to Xq, the properties of the saturated 
turbulence are almost independent of Xq. The saturation level with Xq = 0.4 is slightly higher 
than Xq = 0. The amplitude of time variability is much smaller than those in uniform vertical field 
cases, and the frequency is higher (see Fig. ^). Figure |l6| illustrates the evolution of the magnetic 
energy in model Y4 ((3q = 100, ReMO = 100, and Xq = 0.4) taken along a slice of constant y. 
The contours show the azimuthal component of the magnetic energy with logarithmic spacing, 
and the velocity field is shown by arrows. The growth of disturbances starts from modes with 



large k, typically k x ~ k z ~ 5 for this case (Fig. 1£, top-left panel). At about 10 orbits the MRI 



saturates, and MHD turbulence persists until the end of the simulation at 100 orbits. No growth 
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of the two-channel flow is evident. Small-scale growth of the MRI occurs everywhere and this 
makes small amplitude and frequent time variations. At 20 orbits, the spatial dispersion in the 
magnetic pressure ({SP^g) 1 ^ 2 / (-P ma g) = 1-3) is larger than those in the density and gas pressure 
((5P 2 ) 1/2 /(P) = 0.15 and (5p 2 ) 1/2 /(p) = 0.094). These numbers are almost constant throughout 
the evolution and similar to the zero- net flux B z cases. 

Ohmic dissipation is found to make a large difference in the saturation level of the magnetic 
energy and stress in these pure azimuthal field runs, as it did for models with other field 
configurations. Because the most unstable mode for nonaxisymmetric perturbation has a larger 
k z than the axisymmetric mode on a vertical field, the effect of ohmic dissipation is expected to 



be important for smaller 77, or a larger ReMO- Figure 17 shows the time evolution of the poloidal 
component of the magnetic energy ((B 2 + B 2 )/8ir) for models with different magnetic Reynolds 
number ReMO- The magnetic field strength and the Hall parameter are /5o = 100 and Xq = 0.4 
for all the models in this figure. The initial magnetic field is purely toroidal (B 2 /87t)/Po = 0.01. 
Initially, the poloidal field is generated by small perturbations from the toroidal field, and thus 
much smaller than the azimuthal component ((B 2 + B 2 )/8it)/Pq ~ 10~ 7 . 

When ReMO > 30 the magnetic energy is amplified by the MRI in the linear phase, with the 
linear growth starting earlier with larger ReMO- The saturation level of the poloidal field energy is 
comparable to the initial toroidal field energy. The toroidal field is also amplified by an order of 
magnitude, and it dominates the other components. The Rcmo = 100 and 30 runs evolve to the 
same saturation level, and the Maxwell stress ((wm))/Po is of the order of 0.01. When ReMO = 10, 
the amplified magnetic energy is sustained, but the saturation level of the poloidal field energy 
((B 2 + B 2 )/8ir)/Po is 3 orders of magnitude smaller. The toroidal field is not amplified at all, 
so that the total magnetic energy is unchanged from its initial value. If the magnetic Reynolds 
number is less than 3, the initial perturbations are decaying until the end of the calculation. 

The time- and volume-averaged properties of the turbulence driven by the MRI in the 
azimuthal field runs are listed in Table §. Models Y2 (X = 0) and Y4 (X = 0.4) show little 
effect from ohmic dissipation (ReMO = 100). The more resistive (ReMO = 10) models Y8 (Xq = 0) 
and Y10 (Xq = 0.4) are also listed in the table. Quantities in models with the same ReMO are 
similar, which suggests that the nonlinear effect of the Hall term is small for uniform B y models. 
In fact, the effective Hall parameter is small ((X e ^)) = —0.011 and —0.0012 for models Y4 and 
Y10. For the ReMO = 100 runs, the turbulence is significant and the stress a is larger than 0.01. 
The magnetic energy distribution ({B 2 }} : ((B 2 )) : ((B 2 )) Kt 3 : 22 : 1 and the perturbed kinetic 
energy distribution ((pv 2 )) : ((pSv 2 )) : ((pv 2 )) w2:2:l are similar to those in active turbulence 
of the other initial field geometries. For the ReMO — 10 runs, the poloidal field energy is much 
smaller than the toroidal field energy and the ratio is ((B 2 )) / ((B 2 )) ss 3 x 10 3 . The kinetic energy is 
comparable to the poloidal field energy. The stress a is about 10 -5 so that the angular momentum 
transport is inefficient for these models. 
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5.1. Critical Magnetic Reynolds Number 

Here we consider the critical magnetic Reynolds number for a disk with a toroidal field. 



Figure 18 shows the saturation level of the Maxwell stress as a function of the initial magnetic 
Reynolds number Rcmo- All the models in this figure are without the Hall term, only the effect of 
ohmic dissipation is included. The results with different field strength 0q = 100, 400, and 1600 are 
shown. The critical magnetic Reynolds number for significant stress (e.g., {(wm))/Pq > 0.01) is 
Re-M0,cr\t ~ 30, and this is independent of the initial field strength 0q. The saturation level is also 
independent of 0q, and typically {(wm))/Pq ~ 0.04. For ideal MHD, the saturation level is close to 
this value (Hawley et al. 1995). 

If R&mo is l ess than 30, the evolution of the disk is quite different. Ohmic dissipation 
suppresses the MRI, and the stress in the nonlinear regime is very small and of the order of 10 -6 
- 10~ 5 . The timescale to reach the saturation depends on both the initial field strength and the 
magnetic Reynolds number. When ReMO = 30, for example, the saturation occurs at 25 orbits for 
the fa = 100 run (Y5) but it takes more than 200 orbits for the O = 1600 run (Y24). 

Next we consider the effect of the Hall term on the critical magnetic Reynolds number. Figure 
19 shows the saturation level of the poloidal magnetic energy {((B% + B^)/8tt))/Pq for models 



including the Hall effect. The initial field strength is the same (0q = 100) for all the models. 
Open circles depict the saturation level in models without the Hall term. When Reuo > 30, 
the saturation level is independent of ReMO, and the poloidal field energy is of the order of 0.01. 
If ReMo < 30, the turbulence in the nonlinear regime is reduced by ohmic dissipation, and the 
saturation level drops dramatically. Filled triangles and circles denote the Xq = 0.2 and 0.4 runs, 
respectively. We find that the critical magnetic Reynolds number is not affected by the size of 
the Hall parameter. Even for very large Xq cases (Xq = 10 and 100), the behavior at ReMo < 10 
is unchanged. Therefore, the critical value for uniform B y models is fleMo,crit ~ 30 for any 0q 
and Xq. This is an order of magnitude larger than for uniform B z models. For comparison, the 
magnetic Reynolds number Re' M0 used in Fleming et al. (2000) is also shown in this figure. The 
critical value for Re' M0 is about 3 x 10 3 , however this value should depend on the initial field 
strength. 



6. DISCUSSION 

6.1. Critical Magnetic Reynolds Number at the Nonlinear Stage 

The critical magnetic Reynolds number discussed in sections 3.4, 4.1, and 5.1 is defined 
using the magnetic field strength in the initial state. However, the initial magnetic field in 
accretion disks is highly uncertain, and only the field strength in the nonlinear, saturated state 
of the instability may be observable. We find that the magnetic Reynolds number defined using 
the time- and volume-averaged vertical magnetic field strength in the nonlinear regime, that is 
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((v\ z /rjQ)), characterizes the saturation amplitude of the MRI very well. Tables | - §, §, and @ 



list ((v^/rjQ,)) for all models calculated in this paper. Figure 2C shows the saturation level of the 
stress a = ((w xy ))/Po as a function of ((v\ z /r]Q)} for all the models. Circles, triangles, and squares 
denote models started with a uniform vertical, zero-net flux vertical, and a uniform toroidal field, 
respectively. Filled marks denote models that include the Hall term. 

The dependence of a on the magnetic Reynolds number is clearly evident. A stress larger 
than 0.01 requires ((v\ z /r]U)) > 1. The solid arrows denote the saturation levels obtained for ideal 
MHD (v\/r]VL — > oo), where the upper and lower arrows (a = 0.29 and 0.044) are the average 
of uniform B z runs and uniform B y runs, respectively, taken from Hawley et al. (1995). The 
saturation level when ((v\ z /r]Q)) > 1 is almost the same as the ideal MHD cases, and thus the 
stress is nearly independent of the magnetic Reynolds number in this regime. However, when 
{{v\ z /rjVt)) < 1, the stress decreases as the magnetic Reynolds number decreases. If turbulence 
cannot be sustained by the MRI, the magnetic energy and stress both decrease in time. In some 
models with {(v\ z /r]fl)) < 1, the magnetic energy is decaying at the end of the calculation, and the 
system is evolving toward the direction ((w xy )) oc ((v\ z /rjfl)) shown by dashed arrow in the figure. 

From Figure 20, the saturation level of the stress is approximately given by 



a~aMRimin^l, -^J , (15) 

where «mri is the stress in the ideal MHD limit (v\ z /r]Q S> 1). For a vertical field with zero- net 
flux or a purely toroidal field, the dispersion in the saturation level is quite small: from Tables |6| 
and |8| the averages of the stress are nearly the same ccmki = 0.0395 ± 0.0120 and 0.0372 ± 0.0106 
for these runs. On the other hand, when the disk has net-flux of B z , the saturation levels are less 
uniform because they depend on the vertical field strength, the vertical box size (Hawley et al. 
1995), and the Hall parameter. The stress is ckmri ~ 0.01 — 1, and can be larger than zero-net 
vertical flux models. The stress obtained in Fleming et al. (2000) is shown by crosses in this 
figure; our results are consistent with this earlier work. 

Figure ^ shows that the vertical field is a key quantity for predicting the saturation amplitude 
of the MRI. Nonaxisymmetric modes of the MRI are unstable when azimuthal field is present, 
but because these modes have a large wavenumber in the vertical direction, they are suppressed 
by a smaller resistivity than axisymmetric modes. Therefore the suppression of the axisymmetric 
instability, which depends on the vertical field strength, determines the critical magnetic Reynolds 
number. This is true so long as the azimuthal field gives an Alfven speed that is subthermal 
VA y < c s (Blaes & Balbus 1994), which is always true in our calculations. Thus, the effect of 
nonideal MHD on the saturation amplitude of the MRI is well characterized by a magnetic 
Reynolds number defined as v\ z /t]Q, with a value larger than unity required for turbulence and 
significant transport. 

The critical value ((v\ z /riQ)) ~ 1 is independent of the initial field strength and geometry, 
and also the Hall parameter. The toroidal component of the magnetic energy is dominant, and 
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an order of magnitude larger than the vertical component in the turbulence driven by the MRI. 
The ratio ((B 2 )) / ((B 2 )} ~ 30 when ((v\ z /r]Q)) > 1. Thus the critical value can be written as 
((vj^/rjQ)) ~ 30. If the magnetic Reynolds number vj^/rjQ is larger than about 30, that means the 
MRI would be operating in the disk. 

6.2. Limit on the Parameter Range of Numerical Simulations 

The constraint on the time step in the numerical algorithm used here (see Paper I) is inversely 
proportional to the initial Hall parameter and the field strength. Thus, when the saturation level 
of the magnetic field is high, and when the Hall parameter is large, this constraint is very severe. 
Thus, we have computed only a few models with large Xq and large Rcmo- 

In paper I, we estimated the ratio of the Hall parameter to the magnetic Reynolds number in 
a weakly ionized gas composed by ions, electrons, and neutrals. When the Hall effect is dominant, 
the Hall parameter takes on values Re~^ <J Xq <J 100i?e^Q. We have computed models with 
small i?ejv/o(~ 0.1) and with the maximum value of Xo(~ 1000) which show inefficient angular 
momentum transport in this case. Thus, we conclude that increasing the Hall parameter Xq 
within the allowed range of values will not change the lower limit of the critical ReMo- In fact, 
with large Xq (model S16) the critical wavelength of the MRI is larger than the scale height of the 
disk, so that growth of the MRI cannot be expected for such a situation. 

For the Rej^o > 1 runs, we examined the effect of a Hall parameter that is of the order of 
unity, the maximum value allowed when ReMo ^ 100. However, if 1 <J Repio <^ 10, the Hall 
parameter Xq can take any value between 10 - 100. We have found in this work that a larger 
Hall parameter enhances the saturation level of the stress, especially in the case of a uniform B z . 
But note that even without the Hall effect, significant turbulence is sustained and the stress a is 
more than 0.01 in this regime. Thus, other effects (such as stronger vertical fields) could be as 
important as the Hall effect for enhancing the stress. 

6.3. Application to Protoplanetary Disks 

At the low temperatures expected in protoplanetary disks, thermal ionization is inefficient, 
except for the innermost regions within about r ~ 0.1 AU of the central star. The dominant 
ionization sources in this case are nonthermal processes, such as X-rays, cosmic rays, and 
radioactive elements. By definition, protoplanetary disks contain dust grains that eventually 
will agglomerate into planetesimals. Because recombination process on the surface of grains can 
be important, the number density and size distribution of dust grains has a large effect on the 
ionization fraction of the gas. Unfortunately, there are many uncertainties regarding dust grains 
in protoplanetary disks, and the characteristics of the grains vary in time due to evolutionary 
effects. For example, grains may grow in size through mutual collisions, while the abundance of 
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grains may decrease due to sedimentation toward the midplane. Thus, the effect of dust grains on 
the ionization state of the gas could be reduced in disks in the late stages of evolution. Therefore, 
it is important to investigate the ionization state of the disk in situations both with and without 
dust grains. 

Assuming there are no dust grains in the disk, the distribution of the Hall parameter and 
the magnetic Reynolds number were calculated at the midplane of some disk models in Paper I. 
Inside r# ~ 80 AU, the Hall parameter is \X\ > 1 and the Hall effect is the dominant nonideal 
MHD effect when the field strength is c s /va = 10 in the disk. Ohmic dissipation is more efficient 
in higher density (inner) regions. The magnetic Reynolds number is less than unity at r < ro ~ 6 
AU. The critical radius ro depends on the field strength in the disk, e.g., ro ~ 2 AU and ru ~ 12 
AU for c s /v\ = 1. In the outer disk r > rjj, both the Hall effect and ohmic dissipation are 
unimportant, and thus the disk is unstable to the MRI. 

According to the results of the nonlinear simulations presented in this paper, angular 
momentum transport is inefficient (a <C 0.01) when ReMO <S 1 for any initial field strength (3q 
and for any size of the Hall parameter Xq. Therefore, at the region 0.1 AU <^ r ^ 6 AU, the 
MRI is suppressed by ohmic dissipation and this region probably forms a "dead zone" (Gammie 
1996). The critical magnetic Reynolds number -ReMo.crit for the onset of active turbulence is 1 
- 30 depending on the strength and geometry of the magnetic field. However, because ReM is a 
steeply decreasing function of r, the corresponding uncertainty in the critical radius ro is at most 
a factor of 3 (see Fig. 2 in Paper I). 

In the region ro < r < rn, the Hall parameter is larger than unity, and thus the Hall term 
can affect the evolution of the MRI. Moreover ohmic dissipation is too small to suppress the 
MRI. Although the nonlinear behavior of the MRI depends on the field geometry, the MRI will 
operate for most of the cases in this region. If there is no net flux in the vertical field (zero-net 
flux B z or uniform B y models), MHD turbulence is initiated by the MRI and sustained for the 
values of ReMO and Xq in this region. Although this region is unstable even without the Hall 
term, the stress could be enhanced due to the Hall effect. If the disk is threaded by a uniform 
vertical field oriented in the opposite direction to the angular velocity vector i?, i.e., Xq <C —1, 
the linear growth of the MRI is suppressed by the Hall effect, forming a "dead zone" with no 
angular momentum transport via magnetic stress. Since suppression of the MRI requires the Hall 
parameter is Xo < —4 everywhere, this may be difficult in actual accretion disks. If the vertical 
field is oriented in the same sense as even in only a small part of the disk, the MRI can grow 
from this region. The unstable region could spread wider and eventually fill the entire region, as 
demonstrated in the zero- net flux B z simulations. 

In summary, we expect the outer regions of protoplanetary disks r > ro to be unstable to 
the MRI, with angular momentum transported effectively by Maxwell stress. Inside of ro is a 
dead zone unless the temperature is T ;> 10 3 K. The critical radius ro is a few AU and this is 
determined mainly by the ohmic dissipation. The typical size of protoplanetary disks is about 100 
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AU. Most of the disk is unstable if dust grains are not present due to settling, or are too large to 
affect the ionization fraction. The critical radius ro, which determines where the disk is unstable 
to the MRI and therefore turbulent, depends strongly on the size distribution and abundance of 
dust grains (Sano et al. 2000). The extent to which grains are mixed vertically through the disk in 
turn depends strongly on whether the gas is turbulent. Hence the evolution of gas and dust grains 
must be solved simultaneously in order to make a consistent scenario of grain settling, growth, 
and ultimately planet formation. 

7. SUMMARY 

The Hall effect on the nonlinear evolution of the MRI has been investigated using local 
3D nonideal MHD simulations. Various models with different field strengths (3q, magnetic 
Reynolds number ReMo = v A0 /r]Q, and Hall parameter Xq have been computed. Our findings are 
summarized as follows. 

1. For uniform B z models, a positive (negative) Xq enhances (suppresses) the nonlinear 
turbulence in the disk generated by the MRI. When Xq > 0, the nonlinear evolution 
shows recurrent appearances of the channel flow, and the saturated Maxwell stress is larger 
than that in the negative Xq case. The saturation level of the magnetic energy and stress 
decreases dramatically when Re mo < ReM0,crit ~ 1- This critical value is independent of 
both the field strength and the size of the Hall term. 

2. For zero-net flux B z models, disorganized MHD turbulence is sustained in the nonlinear 
regime without the growth of the channel flow. The critical magnetic Reynolds number 
ReM0,crit for the initial state depends on the initial field strength (3q and the Hall parameter 
Xq, and is in the range of 1 - 30 for our models. The Hall effect can reduce the critical value 
R&M0,crit, however the difference is at most an order of magnitude. 

3. For uniform B y models, the effect of the Hall term can be seen only during the linear growth 
phase of the instability. Disorganized turbulence lasts more than hundred orbits when the 
initial magnetic Reynolds number is ReMo ^ 30, and this critical value is independent of 
both the field strength and the size of the Hall term. The characteristics of the saturated 
turbulence are quite similar to those in the zero-net flux B z models. 

4. The condition for turbulence and significant transport in the nonlinear regime is found to 
be given by ((v Az /rjQ)) ^> 1 (or ((v a /t]Q)) ^ 30), where vaz is the Alfven speed computed 
from only the vertical field component in the nonlinear regime. This is independent of the 
strength and geometry of the initial magnetic field and the Hall parameter. If the magnetic 
field strength in a disk is estimated observationally and the magnetic Reynolds number 
v\/r]Q, is larger than about 30, this would imply the MRI is operating in the disk. 
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We have applied the results of our simulations to the dynamics of protoplanetary disks. We 
conclude the stability of such disks is determined mainly by the distribution of the magnetic 
Reynolds number. The Hall effect makes little change in the critical value of Reu- The size of the 
dead zone where the MRI may be suppressed by ohmic dissipation is sensitive to the characteristics 
of dust grains in the disk. When small dust grains are well-mixed vertically in the disk, the dead 
zone extends to a few tens of AU from the central star. However, the dead zone is considerably 
smaller if small grains are not present. 

Computations were carried out on VPP300/16R and VPP5000 at the National Astronomical 
Observatory of Japan and VPP700 at the Subaru Telescope, NAOJ. This work is supported by a 
grant from the NASA OSS program. 
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Table 1. Uniform B z Simulations (Standard Box Size and Resolution 21 ) 



Model 


0o 




X 


Acrit / H 




Orbits 


((w M ))/Po 






Zl 


3200 


100 


4 


0.14 


0.75 


25 


0.181 


0.0348 


5.92 x 10 3 


Z2 


3200 


100 


2 


0.11 


0.75 


50 


0.162 


0.0338 


4.41 x 10 3 


Z3 


3200 


100 





0.064 


0.74 


50 


0.108 


0.0218 


2.84 x 10 3 


Z4 


3200 


100 


-2 


0.0 


0.70 


50 


0.0394 


0.0124 


1.23 x 10 3 


Z5 


3200 


10 





0.064 


0.70 


25 


0.0896 


0.0227 


271 


Z6 


3200 


1 


4 


0.15 


0.57 


25 


0.184 


0.0268 


53.0 


Z7 


3200 


1 


2 


0.12 


0.51 


25 


0.116 


0.0237 


32.3 


Z8 


3200 


1 





0.091 


0.43 


25 


0.0564 


0.0140 


15.6 


Z9 


3200 


1 


-2 


0.064 


0.28 


25 


0.0272 


0.00956 


7.72 


Z10 


3200 


0.3 





0.22 


0.20 


50 


0.0187 


0.00426 


1.48 


Zll 


3200 


0.1 


4 


0.48 


0.14 


100 


5.05 x 10~ 4 


0.00179 


0.138 


Z12 


3200 


0.1 


2 


0.54 


0.11 


100 


1.16 x 10~ 4 


3.96 x 10~ 4 


0.108 


Z13 


3200 


0.1 





0.64 


0.074 


100 


3.07 x 10~ 4 


0.00136 


0.122 


Z14 


3200 


0.1 


-2 


0.91 


0.037 


200 


6.44 x 10~ 6 


1.74 x 10" 5 


0.101 


Z15 


800 


1 


2 


0.25 


0.51 


25 


0.892 


0.0927 


38.7 


Z16 


800 


1 





0.18 


0.43 


25 


0.178 


0.0322 


9.00 


Z17 


800 


1 


-2 


0.13 


0.28 


25 


0.0323 


0.0117 


3.12 


Z18 


12800 


100 





0.032 


0.74 


25 


0.0251 


0.00734 


2.54 x 10 3 


Z19 


12800 


10 





0.032 


0.70 


25 


0.0340 


0.00944 


375 


Z20 


12800 


1 


2 


0.061 


0.51 


25 


0.0350 


0.00924 


39.4 


Z21 


12800 


1 





0.045 


0.43 


25 


0.0262 


0.00734 


27.9 


Z22 


12800 


1 


-2 


0.032 


0.28 


25 


0.0133 


0.00399 


12.1 


Z23 


12800 


0.3 





0.11 


0.20 


50 


0.0147 


0.00391 


3.86 


Z24 


12800 


0.1 





0.32 


0.074 


100 


4.09 x 10~ 4 


5.89 x 10" 5 


0.127 



a Box size is H x AH x H and grid resolutions is 32 x 128 x 32. 



Table 2. Uniform B z Simulations (High Resolution 3 ") 



Model 


Po 


ReMo 


X() Acrit /H 


^max / ^ 


Orbits 


((w M ))/Po 


({wr))/Po 




Z2H 


3200 


100 


2 0.11 


0.75 


10 


0.128 


0.0262 


5.36 x 10 3 


Z4H 


3200 


100 


-2 0.0 


0.70 


10 


0.0301 


0.0102 


1.35 x 10 3 



a Box size is H x 4H x H and grid resolutions is 64 x 256 x 64. 
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Table 3. Uniform B z Simulations (Large Box Size a ) 



Model 


A> 


ReMo 


Xq A cr it/-ff 


•^rnax / ^ 


Orbits 


«™m))/Po 






Z2L 


3200 


100 


2 0.11 


0.75 


25 


0.159 


0.0428 


6.91 x 10 3 


Z4L 


3200 


100 


-2 0.0 


0.70 


25 


0.0676 


0.0215 


2.47 x 10 3 



a Box size is 2H x 8H x 2H and grid resolutions is 64 x 256 x 64. 



Table 4. Time- and Volume-Averaged Values in Uniform B z Simulations 



Quantity 


Z2 


Z3 


Z4 


Z12 


Z13 


Z14 


Po 


3200 


3200 


3200 


3200 


3200 


3200 


ReMo 


100 


100 


100 


0.1 


0.1 


0.1 


X 


2 





-2 


2 





-2 


{{Bl/^))/P 


0.0435 


0.0272 


0.00974 


3.79 x 10~ 5 


1.16 x 10~ 4 


1.89 x 10~ 6 


K/87r))/Po 
{(B*/8n))/P 


0.314 


0.210 


0.0742 


2.21 x 10~ 4 


6.64 x 10~ 4 


2.22 x 10~ 5 


0.0139 


0.00913 


0.00386 


3.37 x 10~ 4 


3.81 x 10~ 4 


3.14 x 10~ 4 


((pv 2 x /2))/P 


0.0513 


0.0333 


0.0193 


0.00154 


0.00335 


5.09 x 10~ 5 


((pdv 2 y /2))/P 


0.0458 


0.0301 


0.0125 


0.00129 


0.00291 


2.19 x 10~ 5 




0.0160 


0.0135 


0.00816 


0.00833 


0.00812 


0.00514 


((P))/Po 


2.78 


4.63 


2.68 


26.6 


84.7 


1.52 


((wm))/((w b )) 


4.78 


4.95 


3.19 


0.294 


0.226 


0.371 


({w M ))/{{B*/tor)) 


0.454 


0.454 


0.470 


0.457 


0.405 


0.297 


«P»/«B 2 /8tt» 


7.81 


19.5 


32.0 


1.04 x 10 5 


1.12 x 10 5 


7.03 x 10 4 


0.318 


0.324 


0.477 


43.8 


19.0 


241 


{{Bl))/{{Bl)) 


3.12 


2.98 


2.52 


0.112 


0.304 


0.00600 


{{B 2 v ))/{Bl)) 


22.5 


23.0 


19.2 


0.657 


1.74 


0.0707 


((v 2 J/((v 2 )) 


3.21 


2.46 


2.36 


0.184 


0.413 


0.00992 


((/?)) 


2.87 


2.23 


1.54 


0.155 


0.359 


0.00426 


104 


316 


448 


6.78 x 10 4 


1.71 x 10 5 


4.55 x 10 3 


((Re M )) 


1.27 x 10 5 


8.02 x 10 4 


2.88 x 10 4 


0.191 


0.371 


0.108 




0.00277 





-0.0344 


1.41 





-1.86 




0.195 


0.130 


0.0518 


5.12 x 10~ 4 


0.00167 


2.38 x 10~ 5 
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Table 5. Zero-Net Flux B z Simulations (Standard Box Size and Resolution 8 ") 



Model 






^0 


A C rit/-ff 


^"iriax / ^ 


Orbits 


({w M ))/P 


{{wr))/P 






kj 1 


OZiUU 


1 nn 


n 
u 


n nfi/i 


74 


^.n 

uu 


n m 


n ni ^n 




uou 


S9 


OZiUU 


1 nn 


9 


n 1 1 




^.n 

UU 


n n9i ^ 

U. UZi 1 o 


n nnfi7^ 




Ul ( 




OZiUU 


i nn 


A 


n i a 

U. 14: 


7^ 


^.n 

UU 


n n^9n 

U. UUZiU 


n ni w 




1 Q9 y in 3 

i . y Zi a i u 




OZiUU 


i n 

1U 


n 

u 


n nfizt 

u. uu^ 


70 


^.n 

UU 


n m 77 

U.U 1 1 1 


n nns^ 




77 A 




OZiUU 


i n 

1U 


9 
i* 


n 1 1 


79 


^n 

UU 


n Vf>AA 

U. UZi44: 


n nnfi77 




uo . o 


Sfi 




1 

1U 


4 


14 

U. 14: 


7^ 




n n497 


n m 1 8 




1 66 

1UU 


b I 


oonn 
o/UU 


Q 

o 


U 


U.Uoo 


n fin 


Kn 


z.Ul X 1U 


Z.Zo X 1U 


6 


U.UU4U1 


S8 


3200 


1 





0.091 


0.43 


50 


3.95 x 10~ 6 


3.81 x 10" 


6 


2.63 X 10~ 4 


S9 


3200 


1 


2 


0.12 


0.51 


50 


8.83 x 10~ 7 


1.03 x 10" 


6 


6.66 x 10~ 6 


S10 


3200 


1 


4 


0.15 


0.57 


50 


0.0326 


0.00881 




13.4 


Sll 


3200 


0.3 





0.22 


0.20 


50 










S12 


3200 


0.3 


2 


0.21 


0.28 


50 










S13 


3200 


0.3 


4 


0.21 


0.34 


50 










S14 


3200 


0.3 


100 


0.65 


0.70 


50 


7.02 x 10~ 5 


1.86 x 10" 


7 


0.00629 


S15 


3200 


0.1 


100 


0.66 


0.62 


50 










S16 


3200 


0.1 


1000 


2.0 


0.74 


50 










S17 


800 


100 





0.13 


0.74 


50 


0.0173 


0.0106 




194 


S18 


800 


30 





0.13 


0.73 


50 


0.0159 


0.0130 




59.8 


S19 


800 


10 





0.13 


0.70 


50 


1.41 x 10~ 5 


2.44 x 10" 


6 


7.37 x 10~ 4 


S20 


800 


3 





0.14 


0.60 


50 


1.74 x 10~ 7 


3.29 x 10" 


6 


1.75 x 10~ 5 


S21 


12800 


100 





0.032 


0.74 


100 


0.0265 


0.0134 




4.31 x 10 3 


S22 


12800 


10 





0.032 


0.70 


100 


0.0239 


0.0152 




358 


S23 


12800 


3 





0.034 


0.60 


100 


0.0349 


0.0153 




160 


S24 


12800 


1 





0.045 


0.43 


100 


9.90 x 10~ 6 


4.15 x 10" 


6 


0.00179 


S25 


12800 


0.3 





0.11 


0.20 


100 




1.27 x 10" 


6 


1.66 x 10~ 6 



a Box size is H x AH x H and grid resolutions is 32 x 128 x 32. 
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Table 6. Time- and Volume- Averaged Values in Zero- Net Flux B z Simulations 



Quantity 


SI 


S3 


S8 


S10 


Average 3, 


0o 


3200 


3200 


3200 


3200 




ReMo 


100 


100 


1 


1 




X 





4 





4 




((B 2 /8TT))/Pn 
W-^xl u " III 1 U 

((B 2 v /8n))/P 


0.00584 


0.0138 


1.11 x 10~ 8 


0.00938 


0.00821 ± 0.00274 


0.0404 


0.115 


3.71 x 10~ 4 


0.0743 


0.0604 ±0.0237 


((B 2 /8n))/P 


0.00207 


0.00605 


8.24 x 10~ 8 


0.00416 


0.00339 ±0.00129 


((pv 2 x /n/Po 


0.0281 


0.0206 


1.28 x 10~ 4 


0.0146 


0.0266 ±0.0119 


<ip5vl/2))/P a 
((pv 2 z /2))/P 


0.0199 


0.0143 


9.69 x 10~ 7 


0.0102 


0.0176 ±0.0075 


0.00686 


0.00766 


7.95 x 10~ 7 


0.00565 


0.00824 ±0.00304 


((P))/Po 


2.58 


1.94 


1.04 


3.32 




((wm))/((wr)) 


1.29 


3.91 


1.04 


3.70 


2.51 ±0.98 


((w M ))/((By8n)) 


0.362 


0.403 


0.0106 


0.390 


0.394 ±0.018 


((P))/((By8n)) 
{{ P 5v*/2))/{{B*/8k)) 


55.5 


15.0 


2.81 x 10 3 


39.8 




1.18 


0.330 


0.349 


0.364 


0.869 ±0.447 


((B 2 ))/((B 2 )) 


2.83 


2.28 


0.134 


2.26 


2.46 ±0.22 


((B 2 y ))/((B 2 )) 


19.6 


19.0 


4.51 x 10 3 


17.9 


18.0 ±3.1 


»/«0> 


4.10 


2.69 


161 


2.58 


3.15 ±0.46 


(H))/((v 2 )) 
«/?» 


2.90 


1.86 


1.22 


1.80 


2.10 ±0.32 


1.23 x 10 3 


204 


2.81 x 10 5 


562 




((Re M )) 


1.52 x 10 4 


4.51 x 10 4 


1.19 


288 




((X eS )) 





-0.0181 





-0.0335 




a 


0.0298 


0.0653 


7.76 x 10~ 6 


0.0414 


0.0385 ±0.0120 



a Models with ((v% z /riQ.)) > 1 are considered for the average (SI - S6, S10, S17, S18, and S21 
- S23). 
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Table 7. Uniform B y Simulations (Standard Box Size and Resolution* 1 ) 



Model 




ReMo 


X 


AMRi/-ff b 


Orbits 


{(w M ))/Po 






VI 
I 1 


i nn 


1UUU 


n 
U 


n £Q 
U.Oo 


1 nn 
1UU 


n nQ^ft 
U.UoOD 


a nno co 
U.UUoOZ 


Zyo 


i Z 


i nn 


1 nn 
1UU 


(J 


n £;q 
U.Oo 


1 nn 
1UU 


U.Uooo 


n nnooo 
U.UUsoo 


oq ft 
zo.D 


Yo 


i nn 


1 nn 

1UU 


n o 
U.z 


n £:q 
U.Oo 


1 nn 
1UU 


n nQft a 


n nno7n 
U.UU5 (V 


oZ.l 


I 4 


i nn 


1 nn 
1UU 


n a 
U.4 


n £Q 
U.Oo 


1 nn 
1UU 


n n^/i 1 
U.U041 


n m ok 
U.Ulzo 


Oo.Z 


I 


i nn 


on 
oU 


n 
U 


n £Q 
U.Oo 


1 nn 
1UU 


n no^ft 
U.UZoo 


n nn^i £. 
U.UUOiO 


ft 1 A 

0.14 


I 


i nn 


OA 

oU 


n o 
U.Z 
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Y10 


100 


10 


0.4 


0.63 


100 


4.39 x 10~ 5 


1.48 x 10~ 5 


0.00636 


Yll 


100 


10 


10 


0.63 


200 


3.41 x 10~ 7 




6.57 x 10~ 4 


Y12 


100 


3 





0.63 


100 








Y13 


100 


3 


0.2 


0.63 


100 








Y14 


100 


3 


0.4 


0.63 


100 








Y15 


100 


3 


10 


0.63 


100 








Y16 


100 


3 


100 


0.63 


100 








Y17 


400 


1000 





0.31 


200 


0.0279 


0.00606 


816 


Y18 


400 


100 





0.31 


200 


0.0286 


0.00616 


84.6 


Y19 


400 


30 





0.31 


200 


0.0261 


0.00564 


21.4 


Y20 


400 


10 





0.31 


200 


5.34 x 10~ 6 


2.98 x 10~ 6 


5.45 x 10~ 4 


Y21 


400 


3 





0.31 


200 


5.59 x 10~ 7 


1.28 x 10~ 7 


7.83 x 10~ 6 


Y22 


1600 


1000 





0.16 


200 


0.0223 


0.00541 


2.50 x 10 3 


Y23 


1600 


100 





0.16 


200 


0.0189 


0.00467 


200 


Y24 


1600 


30 





0.16 


400 


0.0234 


0.00512 


77.5 


Y25 


1600 


10 





0.16 


200 


3.46 x 10~ 6 


4.51 x 10~ 6 


4.32 x 10~ 4 


Y26 


1600 


3 





0.16 


200 


2.20 x 10~ 6 


4.13 x 10~ 6 


8.51 x 10~ 5 



a Box size is H x AH x H and grid resolutions is 32 x 128 x 32. 

b A M RI = 27TO A o/^ 
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Table 8. Time- and Volume- Averaged Values in Uniform B y Simulations 



Quantity 


Y2 


Y4 


Y8 




Y10 




Average a 


A> 


100 


100 


100 




100 






ReMo 


100 


100 


10 




10 






X 





0.4 







0.4 






II E?2 /q_\\ / r> 


f\ nnonn 

u.uusyu 


U.U141 


o.41 X 1U 


■6 


O.Oo X iU 


-6 


U.UU / 4o ± U.UUz4U 


((B 2 v /Sn))/P 


0.0702 


0.109 


0.00907 




0.00943 




0.0610 ±0.0177 


{{Bl/^))/P 


0.00288 


0.00538 


2.44 x 10- 


■6 


6.32 x 10" 


-6 


0.00247 ±0.00102 


((PV 2 J2))/P 


0.0139 


0.0193 


8.98 x 10- 


■6 


3.10 x 10" 


-5 


0.0117 ± 0.0029 


((pSv 2 y /2))/P 

((pv 2 /n/Po 


0.0105 


0.0157 


2.42 x 10" 


■5 


3.66 x 10" 


■5 


0.00872 ± 0.00259 


0.00583 


0.00833 


1.20 x 10" 


■5 


3.84 x 10" 


-5 


0.00490 ±0.00134 


((P))/Po 


6.44 


5.83 


1.01 




1.02 








4.23 


4.32 


3.74 




2.95 




4.33 ±0.20 


({w m ))/((B 2 /8k)) 
((P))/({B 2 /8n)) 


0.447 
81.7 


0.433 
46.6 


0.00245 
111 




0.00465 
108 




0.442 ±0.022 


{{p5v 2 /2))/{{B^/^)) 


0.383 


0.347 


0.00498 




0.0112 




0.372 ±0.018 


((B 2 J/((B 2 )) 


3.09 


2.61 


1.40 




0.891 




3.10 ±0.28 


((B 2 ))/((B 2 )) 


24.4 


20.3 


3.73 x 10 3 




1.49 x 10 3 




25.6 ±2.5 


((v 2 J/((v 2 z )) 


2.38 


2.31 


0.747 




0.806 




2.42 ±0.13 


({Sv 2 J/{{v 2 )) 


1.80 


1.89 


2.02 




0.953 




1.78 ±0.10 


1.39 x 10 3 


688 


114 




115 






((Re M )) 


825 


1.29 x 10 3 


9.08 




9.45 






((X**)) 





-0.0112 







-0.00119 






a 


0.0436 


0.0666 


2.82 x 10" 


■5 


5.88 x 10" 


-5 


0.0372 ±0.0106 



a Models with ((v% z /r)Q,)) > 1 are considered for the average (Yl - Y7, Y17 - Y19, and Y22 - 
Y24). 
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Fig. 1. — Time evolution of the volume-averaged magnetic energy {B 2 /87r)/Po for models Z2 
(Xq = 2), Z3 {Xq = 0), and Z4 (Xq = —2). The plasma beta and the magnetic Reynolds number 
of these models are (3q = 3200 and Rcmo = 100. 
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Fig. 2. — Time evolution of the volume-averaged Maxwell stress (—B x B y /4:Tr)/Po for models Z2 
(Xq = 2), Z3 (Xq = 0), and Z4 (Xq = —2). The plasma beta and the magnetic Reynolds number 
of these models are Pq = 3200 and ReMo = 100. 
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Fig. 3. — Slices in the x-z plane at y = —2H and in the y-z plane at x = —0.5H of the magnetic 
energy, log(B 2 /87rPo) (colors), and perturbed velocity field, 5v (arrows), in model Z2 (/3q = 3200, 
Re M0 = 100, and X = 2) at 18, 21, and 25 orbits. 



Fig. 4. — Slices in the x-z plane at y = —2H and in the y-z plane at x = —0.5H of the magnetic 
energy, log(B 2 /87rPo) (colors), and perturbed velocity field, 5v (arrows), in model Z4 (/3q = 3200, 
Re M0 = 100, and X = -2) at 10, 17, and 25 orbits. 
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Fig. 5. — Saturation level of the Maxwell stress as a function of the Hall parameter Xq for the 
models with (5q = 800, 3200, and 12800. The magnetic Reynolds number is Reuo = 1 for all the 
models. 




Fig. 6. — Saturation level of the Maxwell stress as a function of the Hall parameter Xq for the 
models with Rcmq = 100, 1, and 0.1. The plasma beta is (3q = 3200 for all the models. 
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Fig. 7. — Power spectra of magnetic field in models (a) Z2 and (b) Z4 along the k x , k y , and k z 
axes. The wavenumbers are normalized in terms of H/2tt and energies in terms of initial pressure 
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Fig. 8. — Power spectra of (a) magnetic field and (b) velocity in model Z13 along the k x , k y , and k z 
axes. The wavenumbers are normalized in terms of H/2ir and energies in terms of initial pressure 
Pq. The arrows lie at the diffusion scale k<as = vao/v an< l at the critical wavenumber for the MRI 
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Fig. 9. — Saturation level of the Maxwell stress as a function of the magnetic Reynolds number 
ReMO- Open circles and triangles denotes the models without Hall term (Xq = 0) for (3q = 3200 and 
(3q = 12800, respectively. The models including the Hall term are shown by filled circles (Xq = 4) 
and triangles (Xq = — 2). 
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Fig. 10. — Time evolution of the volume-averaged Maxwell stress (— B x B y / An) / Pq for models SI 
(Xq = 0), S2 (Xq = 2), and S3 (Xq = 4). The plasma beta and the magnetic Reynolds number of 
these models are 0o = 3200 and Reuo = 100. 



Fig. 11. — Slices in the x-z plane at y = of the azimuthal component of the magnetic 
energy, \og{B 2 y / '8tt Pq) (colors), and magnetic field vector, B (arrows), in model S3 (/3o = 3200, 
ReMO = 100, and X = 4) at 35, 40, 45, and 50 orbits. 
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Fig. 12. — Time evolution of the volume-averaged magnetic energy {B 2 /8ir) /P$ for models S8 
(Xq = 0), S9 (Xq = 2), and S10 (Xq = 4). The plasma beta and the magnetic Reynolds number of 
these models are (3q = 3200 and Re mo = 1- 
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Fig. 13. — Saturation level of the Maxwell stress as a function of the magnetic Reynolds number 
ReMO for zero-net flux B z models. All the models are without the Hall effect (Xq = 0), and the 
initial field strength is f3 = 800, 3200, and 12800. 
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Fig. 14. — Saturation level of the magnetic energy as a function of the magnetic Reynolds number 
ReMo for zero- net flux B z models (/3q = 3200). Open circles denote the models with only the ohmic 
dissipation (Xq = 0), and the other marks are including also the Hall effect (Xq = 2, 4, 100, and 
1000). 
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Fig. 15. — Time evolution of the volume-averaged Maxwell stress (—B x B y / An) / P$ for models Y2 
(Xq = 0), Y3 (Xq = 0.2), and Y4 (Xq = 0.4). The plasma beta and the magnetic Reynolds number 
of these models are (5q = 100 and Reuo = 100. 



Fig. 16. — Slices in the x-z plane at y = of the azimuthal component of the magnetic 
energy, log(-B^/87rPo) (colors), and perturbed velocity field Sv (arrows) in model Y4 ((3q = 100, 
Re M o = 100, and X = 0.4) at 5, 10, 15, and 20 orbits. 
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Fig. 17. — Time evolution of the volume-averaged poloidal component of the magnetic energy 
{{B 2 X + B*)/8ir)/P for models Y4 (Re M o = 100), Y7 (Re M o = 30), Y10 (Re M0 = 10), and Y14 
(ReMo = 3). The plasma beta and the Hall parameter of these models are (3$ = 100 and Xq = 0.4. 
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Fig. 18. — Saturation level of the Maxwell stress as a function of the magnetic Reynolds number 
Re mo for uniform B y models. All the models are without the Hall effect (Xq = 0), and the initial 
field strength is fa = 100, 400, and 1600. 
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Fig. 19. — Saturation level of the poloidal component of the magnetic energy as a function of the 
magnetic Reynolds number ReMo for uniform B y models (/3b = 100). Open circles denote the 
models with only the ohmic dissipation (Xq = 0), and the other marks are including also the Hall 
effect (X = 0.2, 0.4, 10, and 100). 



-43- 




-6 -4 -2 2 4 6 
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Fig. 20. — Saturation level of the stress ((w xy ))/Po as a function of the magnetic Reynolds number 
{{v\ z /i]V!)) at the nonlinear stage. Circles, triangles, and squares denote the models started with a 
uniform vertical field, zero-net flux vertical field, and a uniform toroidal field, respectively. Filled 
marks means the results of the models including the Hall term, and open marks are without the 
Hall effect. Crosses are the results obtained by Fleming et al. (2000). Upper and lower solid arrows 
indicate the averages of the ideal MHD runs in Hawley et al. (1995) for initially uniform B z and 
initially uniform B y , respectively. 
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